Increased Brownian force noise from molecular impacts in a constrained volume 
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We report on residual gas damping of the motion of a macroscopic test mass enclosed in a nearby 
housing in the molecular flow regime. The damping coefficient, and thus the associated thermal force 
noise, is found to increase significantly when the distance between test mass and surrounding walls 
is smaller than the test mass itself. The effect has been investigated with two torsion pendulums 
of different geometry and has been modelled in a numerical simulation whose predictions are in 
good agreement with the measurements. Relevant to a wide variety of small-force experiments, the 
residual-gas force noise power for the test masses in the LISA gravitational wave observatory is 
roughly a factor 15 larger than in an infinite gas volume, though still compatible with the target 
acceleration noise of 3fms~^ Hz~^" at the foreseen pressure below lO^'' Pa. 

PACS numbers: 95.55.Ym, 05.40. Jc, 07.10.Pz 
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The fluctuation-dissipation theorem [T] states that any 
system with dissipation exhibits fluctuations analogous 
to Brownian motion that can be modelled with an exter- 
nal driving force with power spectral density 



Sf{uj) = 4A:TRe 



dF 
dv 



AkTRe{Z{uj)) (1) 



where k is Boltzmann's constant, T is the temperature of 
the system, and Z{lij) is the mechanical impedance of the 
system. Those equilibrium fluctuations are a fundamen- 
tal limit of precision metrology experiments employing 
macroscopic test-masses nominally in perfect free-fall to 
define geodesic reference frames. 

An unavoidable source of dissipation comes from resid- 
ual gas in the experimental volume, which yields a vis- 
cous, frequency independent, impedance Z{uj) = f3. The- 
oretical calculations for small force experiments [3S] and 
direct observations |5l|6j show the gas damping coefficient 
/3 to be proportional to residual gas pressure P and to 
some effective surface area, which depends on the test 
mass geometry. 

A recent calculation T for a cubic test mass has found 
that the translational damping coefficient is 
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where s is the side length of the test mass, and m is the 
mass of the gas molecules. For rotation, relevant for the 
experiments we will present shortly, the coefficient is 
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This calculation assumes that gas and test mass are in 
thermodynamic equilibrium, that gas molecule collisions 
with the test mass are completely inelastic, with prompt 



stochastic re-emission from the surface with a Maxwell- 
Boltzmann velocity distribution and cosine law angular 
distribution. Importantly, this and other calculations as- 
sume that the test mass is surrounded by an infinite vol- 
ume of coUisionless gas, hence the superscript 00, such 
that a molecule emitted from the test mass surface dis- 
appears into the surrounding volume, and that the mo- 
mentum it imparts on the test mass is uncorrelated with 
any subsequent collisions. This assumption breaks down 
when the distance between the test mass and the sur- 
rounding enclosure is only a fraction of the test mass 
size, as in the geometry studied here, where the gap is 
roughly one tenth of the test mass size. 

Our experimental setup probes small forces relevant 
to the quality of free fall achievable for the test masses 
of the Laser Interferometer Space Antenna (LISA) 
gravitational-wave observatory [51 [S] and its precursor 
LISA Pathfinder [ID], using a torsion pendulum. Exter- 
nal force disturbances are suppressed to a level where the 
thermal noise plays an important role in limiting the force 
sensitivity. Measurements of the torsion pendulum free- 
motion damping show a residual gas contribution several 
times larger than predicted by the infinite- volume model, 
accompanied by an observed noise increase consistent 
with mechanical thermal noise. We present here simple 
arguments, based on the impedance to stochastic molec- 
ular flow through the small gaps around the test mass, 
supported by numerical simulation, that quantitatively 
explain this observed increased dissipation. Enhanced 
gas damping has also been observed in micro-mechanical 
resonators |11] , though the interpretations have employed 
an elastic model of the molecule- wall interaction [12] . 

Two different torsion pendulum geometries are used 
for the experimental investigation. Both suspend a hol- 
low replica of the LISA Pathfinder cubic test mass with 
side length s = 46 mm inside a prototype of the LISA 
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FIG. 1: Sketch of the ITM and 4TM pendulum inertial mem- 
ber surrounded by the GRS. Test mass side is 46 mm. Gaps 
between test mass and GRS are 4.0, 2.9, and 3.5 mm on re- 
spectively the X, y, and z faces. The S2 gaps are 8.0, 6.0, and 
8.0 mm gaps on x, y, and z faces, respectively. 



Pathfinder Gravity Reference Sensor (GRS) [T3] that 
completely surrounds the test mass. All GRS, test mass, 
and other inertial member surfaces are gold coated. In 
the ITM torsion pendulum [T3[[TS] the rotation axis coin- 
cides with the test mass z axis of symmetry, as shown on 
the left of Fig. [ij In the 4TM torsion pendulum [16l [H] 
the test mass enclosed by the GRS is one of four iden- 
tical test masses in a cross shaped configuration, sus- 
pended off-center with respect to the fiber axis with a 
r = 0.1065 m arm, as shown on the right of Fig. [I] The 
opposite test mass is enclosed in a similar capacitive sen- 
sor (S2), featuring larger gaps. The 4TM pendulum em- 
ploys a 50 /im diameter tungsten fiber, while the ITM 
pendulum uses a 40 /im fused silica fiber with much lower 
internal dissipation. Both torsion pendulums facilities 
operate at the temperature of 293 K, controlled within 
O.IK. 

For both pendulums we have measured the torsional 
damping coefficient in the pressure range of roughly 
0.003-3 mPa (mean free paths of order 10 km down to 
10 m), by direct measurement of the amplitude decay 
time T of the free pendulum motion 

where X is the torsion pendulum moment of inertia. The 
decay time is computed as the average of the ratio be- 
tween oscillation amplitude, with each amplitude point 
averaged over three periods, and its time derivative, as es- 
timated by the difference between successive points. Un- 
certainty is computed on statistical basis. 

Measurements were typically performed over 1-3 days, 
using large initial amplitudes of order mrad. No ampli- 
tude dependence of r was observed. The ITM and 4TM 
oscillation periods were roughly 460 and 1400 s, respec- 
tively, and their ring-down times in the limit of zero pres- 
sure were roughly 6 years {Q ~ 1.2 x 10^) and 16 days 
(Q ~ 3.2x 10^). In the case of the 4TM torsion pendulum 
we performed the measurement both with and without 
the GRS around the test mass. 



The pressure in the torsion pendulum vacuum cham- 
bers is measured with ion gauges cross-calibrated with a 
certified pressure gauge whose accuracy is within 15%, 
employing the calibration factor recommended for air. 
While we have not as yet identified the dominant gas 
species in our chamber, measurements with the same 
apparatus have confirmed the pressure-dependent radio- 
metric effect model [l5l [T^ at the 10% level, indicating 
similar systematic uncertainty in the pressure measure- 
ments here. 

The measured damping coefficients, plotted in Figs. [2] 
and |3] are found to depend linearly on pressure. The 
residual zero-pressure intercept, from other dissipation 
mechanisms, has been removed from both datasets, and 
the pressure dependence is summarized in Table |Ij The 
ITM data can be directly compared with the infinite vol- 
ume prediction in Eqn. [3] The 4TM prediction is calcu- 
lated from the rotational and translational contributions, 
/3 = 4 {r'^13^ + Pmt)- The infinite-volume model predic- 
tions are several times smaller than the measured values 
and would foresee no difference caused by the removal of 
the GRS, while here it lowers the total 4TM damping by 
more than a factor two. 





measurement oo model simulation 




dl3/dP\m''s\ dl3/dP\m^s\ 9/3/aP [m» s] 


I. 


(4.8 ± 0.1) X 10"* 1.74 xl0~** 4.24x10"® 


II. 


(9.4±0.2) X lO"'' 1.51x10"^ 7.87 xlO"'' 


III. 


(3.7 ± 0.1) X 10"'' 1.51 X 10"'' 3.27 x 10"'' 


IV. 


(5.7 ± 0.3) X 10"'' 4.60x10"'^ 



TABLE I: Comparison of P pressure dependence, obtained 
from measurement, infinite- volume model, and simulation, in 
different configurations: I. ITM, II. 4TM with GRS, III. 4TM 
without GRS, IV. 4TM GRS contribution, obtained as differ- 
ence between configurations II and III. Simulation uncertainty 
are on the fourth significative digit. 



A series of torque noise measurements, performed with 
the ITM pendulum, has confirmed the Brownian nature 
of the gas damping noise. Measurements were performed 
at six different pressures, from Pq = 0.004 to 15mPa. 
Torque noise was estimated with a dual-readout cross- 
correlation technique [HI [TS]. For each measurement, 
lasting 1-3 days, noise spectra were averaged over 8 to 
20 windows of 25000 s each, with 50% overlap and over 
logarithmically-spaced frequency bins, with uncertain- 
ties estimated as the standard deviation divided by the 
square root of the number of points. To identify the gas 
damping noise, the Pq spectrum was subtracted from the 
five higher pressure measurements, yielding a differential 
spectrum ASn for each A^ = (3{P) - l3{Po). Two such 
spectra are shown in the inset of Fig.|4J The gas damping 
noise is best fit by a white noise model between 0.5 and 
12mHz, where the gas damping noise is resolved. The 
averages of each differential spectrum in this bandwidth 
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FIG. 2: Measured ITM gas damping coefficients /3rot as func- 
tion of residual gas pressure P. Infinite-volume model and 
simulation predictions are also shown. 



are plotted as a function of A/3 in Fig. |4] with a linear fit 
yielding ASn = (1-08 ± 0.03) xAkTA/S. The gas damp- 
ing noise is thus white and within 10% of the Brownian 
noise prediction. 

The excess gas damping and associated force noise, 
well above that foreseen by the infinite-volume model, is 
caused by the nearby GRS walls, as can be understood 
by a simple picture of molecular gas flow. Test mass 
motion along the x direction, with velocity d, requires 
gas flow Q, expressed in Pam^s~i^, from the volume on 
one X face to the other, to restore pressure equilibrium. 
The finite molecular flow conductance C of the channels 
provided by the electrode-housing y and z faces gaps re- 
quires a pressure drop AP to maintain this gas flow. In 
steady state conditions, Q = Ps'^v and thus the pressure 
drop must be AP = Q/C = Ps^v/C. The test mass is 
thus subject to a force, F = APs^, proportional to the 
velocity v, and thus to an increase in the viscous damping 
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Gas flow in a finite conductance is a dissipative process. 

In an equivalent picture from the standpoint of the as- 
sociated force noise, the flow impedance of the channels is 
related to the time needed for molecules to random walk 
from one side of the test mass to the other. This ran- 
dom walk requires repeated impacts of a single molecule 
on the same test mass face and thus introduces correla- 
tions between subsequent collisions, with repeated force 
impulses of the same sign slowing the averaging out of 
the net force and thus increasing force noise. 

To verify our hypotheses, we developed, on the same 
assumptions as the infinite-volume model, a numerical 
simulation of the gas dynamics within a simplified geom- 
etry of a cubic test mass inside a cubic housing. The 
gas is supposed to be composed of a single specimen 
with molecular mass of 30 amu. The simulation is com- 



FIG. 3: Measured 4TM gas damping coefficients /?rot as func- 
tion of residual gas pressure P, with and without the GRS 
surrounding one of the test masses. Infinite-volume model 
and simulation predictions are also shown. 
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FIG. 4: Measured averaged gas damping noise difference ASjv 
in the band 0.5-12 mHz as a function of differential damping 
coefficient A/3, compared with the Brownian noise prediction. 
Inset shows two noise spectra, compared against the predicted 
noise level. 



posed of N steps, simulating each one the dynamics of 
a single molecule. The simulation allows the system to 
evolve for a time At chosen long enough to allow time 
for the molecules to stochastically move many times from 
one side of the test mass to the other. Momentum ex- 
change between test mass and gas molecules is recorded 
at each collision. The gas-damping coefficient /3 is com- 
puted from the variance a\ of each component of the 
total force acting on the test mass: 
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An equivalent approach is used to calculate /3^q° from the 
mean square fluctuation of the torque. 

The simulation has been performed as a function of gap 
size, with the results shown in Fig. [Sj In addition to con- 
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FIG. 5: Gas damping /J'''™ obtained from the numerical sim- 
ulation for different test mass side lengths s and gap sizes d, 
normalized to the infinite- volume model prediction /?°°. 



firming the analytic model in Eqns.[2]and[3]for large gaps, 
the data show the expected increased damping for d < s, 
which approaches a power law (d/s)~'^ for vanishing gap. 
This can be roughly understood with the flow impedance 
arguments presented earlier. For s ^ d, the channel 
conductance is of order C ~ d'^ In (s/d) (kT/m) ' [TO] . 
Substituting into Eqn. [5] and comparing with Eqn. [2] 
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This asymptotically approaches a log slope of —2 for van- 
ishing gap, as observed in the simulation data. 

Fig. [5] also shows simulation results obtained simulat- 
ing the slightly asymmetric gaps of the real GRS geom- 
etry. Finally, we also simulate the 6 mm diameter holes 
present on the x faces, which act as vents that reduce the 
damping pressure head, with a roughly 20% reduction of 
the gas damping obtained in the simulation. 

The measured ITM gas damping coefficient can be di- 
rectly compared with the simulation results, adding in 
small contributions due to minor features of the iner- 
tial member computed with the infinite-volume model. 
The gas damping coefficient for the 4TM torsion pen- 
dulum without the GRS is computed from simula- 
tion results for the S2 geometry and from the infinite- 
volume model for the other three test masses /3 = 
3 (r2/3t7 + /3^^) + r^PI^ + /3,"™. The GRS contribution 
in the 4TM gas damping coefficient is computed from 
simulation results and infinite- volume model predictions 
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Results are sum- 



marized in Table IT] The agreement with measurements 
is within 20%. This is good in respect of the approxi- 
mations in the geometrical model of the GRS contained 
in the simulation. We note that in the 4TM pendulum 
the gas damping coefficient is dominated by the transla- 
tional part, the rotational one accounting for roughly 1% 
of the total, and thus the 4TM measurement validates 
the translational damping relevant to force noise. 



Converting the simulated residual-gas Brownian force- 
noise into acceleration noise for LISA and LISA 
Pathfinder, with a 2 kg test mass, yields 



S}/^ = 1.3 X 10^1^ 
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Eqn. |8] is compatible with the target acceleration noise 
of LISA Pathfinder and LISA, the latter of which will 
operate at pressure below 10^^ p^^ reached by venting 
the GRS to space ^. Gas damping is likely to be the 
dominant pressure-related force noise for LISA. We note 
that, unlike many force noise sources, Brownian noise 
from residual gas does not improve at higher frequencies, 
this is important for future gravitational wave missions, 
such as DECIGO and BBO, whose ambitious sensitivity 
goals near 0.1 Hz and shorter intcrferometry arms require 
even lower acceleration noise limits [20, 21_. Finally, the 
gas-damping question effectively confirms the need, also 
based on electrostatic concerns jl3j, for gaps of at least 
several mm: for a similar sized test mass and gaps of 
0.3 mm, roughly the gap employed in several space ac- 
celerometry missions [55J [35] , the resulting acceleration 
noise from this source alone would be an order of magni- 
tude above the LISA goal. 
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